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Abstract: We consider a lattice gas with general short range interactions and a 
Kac potential J 7 (r) of range -f~ , 7 > 0, evolving via particles hopping to nearest 
neighbor empty sites with rates which satisfy detailed balance with respect to the 
equilibrium measure. Scaling space like 7 _1 and time like 7~ 2 , we prove that in the 
limit 7 — > the macroscopic density profile p(r, t) satisfies the equation 



.(p)v- 



(*) 



Sp(r) 

Here a s (p) is the mobility of the reference system, the one with J = 0, and ^F(p) = 
J [fs(p( r )) ~ \p( r ) J J{ r ~ r')p(r')drdr'}, where f s (p) is the (strictly convex) free 
energy density of the reference system. Beside a regularity condition on J, the only 
requirement for this result is that the reference system satisfy the hypotheses of the 
Varadhan— Yau Theorem [VY] leading to (*) for J = 0. Therefore (*) holds also 
if T achieves its minimum on non constant density profiles and this includes the 
cases in which phase segregation occurs. Using the same techniques we also derive 
hydrodynamic equations for the densities of a two component A-B mixture with 
long range repulsive interactions between A and B particles. The equations for the 
densities pa and pb are of the form (*). They describe, at low temperatures, the 
demixing transition in which segregation takes place via vacancies, i.e. jumps to 
empty sites. In the limit of very few vacancies the problem becomes similar to phase 
segregation in a continuum system in the so called incompressible limit [EP], [DG]. 
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The state of a (one component) macroscopic system in equilibrium can be characterized 
by two numbers, the temperature T(= and the chemical potential A. When T and 
A correspond to a single phase (i.e. there is a unique Gibbs measure) then the particle 
density p(T, A) is constant, i.e. spatially uniform. Given a microscopic dynamics for which 
this Gibbs measure is attractive, a disturbance in this uniform density corresponding to 
a profile po(r), r the space coordinate, at some time to is expected to relax towards the 
uniform density p. In certain types of systems (when the variations in temperature and 
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hydro dynamical flows can be neglected, see below) the relaxation of the density profile 
p(t, r) will occur via the diffusion equation 

|£ = V-[£>Vp], t>t , (1.1) 

where the bulk diffusion constant D = D(p) will generally depend also on the temperature 
T, assumed constant and therefore omitted in the notation. 

Equation (1.1) is a continuity equation for p corresponding to a mass current given by 
Fick's law 

j = -D{p)Vp. (1.2) 

To derive (1.1) from microscopic models it is most convenient to write Fick's law in its 
Onsagerian form, 

3 = -*VA, (1.3) 

where a is the conductivity, or mobility, and X(p) is the local chemical potential; A and 
p are related as in the uniform equilibrium system, i.e. we are in a situation of local 
equilibrium. Comparing (1.2) and (1.3) gives 

D = o/ X (1-4), 

where 

is the compressibility (p being the pressure). In (1.5) we have expressed the chemical 
potential at equilibrium as a function of the density p. The relation (1.4) is sometimes 
referred to as the Einstein relation who first used it to relate the diffusion constant of a 
Brownian particle to its steady state mobility in an external field [E, Sp]. 

Mathematically rigorous derivations of (1.2) - (1.5), with cr given by a Green-Kubo 
formula, have been achieved, via the use of the hydro dynamical (diffusive) scaling limit for 
a variety of microscopic models with fixed short (microscopic) range interactions evolving 
via stochastic dynamics (in which particle number is the only conserved quantity), see 
[KL]W. In all these cases the temperature T is in the uniqueness region of the phase 
diagram, i.e. there is a unique phase for all values of A (or p). 

The situation becomes much trickier when we consider temperatures where there is, for 
some value of A, a coexistence of phases with two (or more) different densities, p\ and 
P2, corresponding to liquid and vapor or to fluid and solid. In such cases the macroscopic 
equilibrium system, with a fixed total number of particles corresponding to an average 
density in the interval (pi, P2), will not have a uniform density. Instead it will be segregated 
into macroscopic regions of density p\ and p2 with shapes minimizing the free energy of 

W There also exists a derivation for one system with a Hamiltonian evolution, namely a 
non-interacting gas of point particles moving among fixed periodic scatterers - the Sinai 
billiard with finite horizons [BS], [ Sp]. 
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the surface between them. An equation of form (1.1) is clearly not appropriate now, in 
fact, for any density p G (pi,pz), x(p) i n (1-4) will be infinite and therefore D is formally 
zero (unless a is also infinite, which can be proven in some cases not to happen). Results 
in these directions have been proven for very particular systems, see the remarkable result 
in [R] where degenerate diffusion in the coexistence region is proven, but the general case 
of the evolution of phase domains for systems with interactions is a real challenge for the 
moment. 

To get around this obstacle and derive macroscopic equations for a system undergoing 
phase segregation, some authors studied the time evolution of the macroscopic density pro- 
file in particle systems interacting via long range (compared to the interparticle spacing) 
Kac potentials, [LOP], [GL1]. The microscopic model considered in [GL1] is a lattice gas 
evolving under a particles hopping (Kawasaki exchange) dynamics which satisfies detailed 
balance (is reversible) with respect to the Gibbs canonical (fixed particle number) equi- 
librium measure with Hamiltonian H at temperature j3~ . H consisted of a sum of two 
terms, a short range part H s which may be thought of as a nearest neighbor interaction 
and a Kac potential H 1 characterized by a range parameter 7 , namely 

#7 = -\ E i dj ln\ x - y\)v(x)v(y), (i-6) 

x,y 

where x,y G A 7 C 7L d (say a torus of diameter [7 -1 ]) and rj(x) = 0,1 specifies the oc- 
cupancy of site x. They argued that in the diffusive hydro dynamical scaling limit, corre- 
sponding to 7 — > with space and time scaling like 7 _1 and 7 -2 respectively, the density 
profile p(r, t) would satisfy a parabolic integro-differential equation of the form, 



dp 
dt 



a s (p)V 



dp 



= V-[a s {p)VX{p)\ 



(1.7) 



where 



T = J f s (p(r))dr—-J J J(r — r')p(r)p(r')drdr' 



(1.8) 



f s is the free energy density and a s is the mobility of the system with only short range 
interactions. 

Equations (1.5) and (1.7) were proven in [GL1] for the case where the short range 
interactions H s consisted only of the hard core exclusion, i.e. no more than one particle 
per site, in which case f s and a s reduce to f® and a®, with 



and 



/ s °(p) = /?[plogp+(l-p) log(l-p)], 



a° s (p) = pp(l-p). 



(1.9) 



(1.10) 



The validity of (1.7) and (1.8) for the case of nontrivial short range interactions was 
conjectured in [§3,GL1] for values of above the critical temperature of the reference 
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(short range) system, when f s (p) is a strictly convex function of p. They gave heuristic 
arguments for the validity of (1.7) with A the local chemical potential of the system with 
the Kac potential, but with the mobility a being the same as what it would be in the 
system without the slowly varying Kac potential. This generalized a similar conjecture by 
Spohn [Sp] for the case of an external long range interaction. Besides the case in [GL1], 
the conjecture was shown to be valid also for some other special cases [BL] , [MM] , [AX] . In 
this paper we prove these conjectures for the general short range interactions in the case 
of Ising spins, that is in the case of a system with Hamiltonian H 1 = H s + H 1 . In fact 
the methods we use here extend in a straightforward manner to general systems for which 
a diffusion equation can be proven in the absence of H 1 . We illustrate this by deriving in 
Section 5 integro-differential equations for a binary mixture which may undergo a demixing 
transition. Our results extend to the case of a system on which a weak external force is 
acting, characterized by the Hamiltonian 

X 

and also to the case of a weak (of order 7) external force which is not the gradient of a 
potential, like a costant force in a torus. 

We note that solutions to (1.7) corresponding to interface dynamics are considered in 
[GL2] . There is also a review of various results about these models, including the case in 
which there is no conservation law (Glauber dynamics) [GLP]. 

Informal Description of Model. To be concrete, we restrict ourselves to the case H s (rj) = 
K J2x y\x- y \=i v( x )v(y)- The Kawasaki dynamics is defined in terms of Poisson jump rates 
depending on the energy differences. An example of such rates is 

cJ XiV) ( V ) = exp (-(3 [H^rf") - £T fa)] /2) , (1.11) 

where x, y are nearest neighbor sites and r/ x ' y is the configuration in which the sites x and 
y exchange their occupation numbers. The microscopic current w Xjy through the bond 
(x, y) is the rate at which a particle jumps from x to y minus the rate at which a particle 
jumps from y to x, namely 

u>2, v (v) = v(x)cJ XiV) (ri) - v(y)c( Xty) (v)\ ■ (1-12) 

Equation (1.12) will determine the form of the macroscopic current in the hydro dynamic 
equation. 

A key ingredient in our analysis is that the dynamics with J 7^ is a weak perturbation 
of the J = dynamics. This can be seen both at the level of the rates and of the current. 
In particular for the current some straightforward expansions with respect to the small 
parameter 7 (see Section 4) lead to 

w 2,x+ ei (v) =c x , x+ei (V) lv(x) ~ V(x + et)] 

-id ( L13 ) 
+ -y <* +ei [v(x) ~ V(x + e,)] 2 (d t J * rj) (x) + 0( 7 2 ), 
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where is the unit vector on the lattice in the direction i, * denotes spatial discrete 
convolution and the superscript on c and w denotes the case in which there is no long 
range force (J = 0). Equation (1.13) indicates clearly that the dynamics with J ^ 
is a weak perturbation of the reference (J = 0) dynamics. The case studied in [GL1] 
corresponds to c° = 1 in (1.13). The local equilibrium expectation of x+e . gives the 
macroscopic flux j for this simplified model in the form of a term due to the exclusion 
dynamics and one due to the mean field force: 

j = - Vp + a(p) V J * p, (1.14) 

where * denotes the (standard) spatial convolution, J*p(r) = J J(r — r')p(r')dr' , and the 

mobility a(p = (3p(l — p) is just the expectation of ((3/2) [rf(x) — rj(x + e^)] . 

The naive extension of the above argument to the general case, i.e. just applying local 
equilibrium to (1.13), would give a macroscopic current of the form 

j = —DVp + a(p) V J * p, (1.15) 

where D is, as before, the diffusion coefficient found in [VY], but a would not satisfy 
the Einstein relation (1.4). This apparent contradiction is explained by the fact that the 
microscopic current of the reference system in the general case, i.e. the first term in the 
right-hand side of (1.13), is in not a lattice gradient of a function of the configuration, i.e. 
there is no function h such that (rj) = h(T y rj) — h(r x rj), where r x is the translation in 
configuration space, unless K = (i.e. H s = 0). The non gradient nature of the dynamics 
is responsible for the presence in (1.15) of a third term which, combined with the second 
term, gives the correct expression for the macroscopic current. 

The main ingredient in our derivation is the recent work of Varadhan and Yau [VY] 
who proved (1.1) for a lattice gas with general short range interactions at small f3. This is 
a major achievement since previous derivations all required that the dynamics be either of 
gradient type, or that the invariant measure be of product type, i.e. independent occupa- 
tion values at different sites, see [KL] for a complete treatment and extended bibliography 
on this topic. The class of non-gradient models with invariant measures of product type 
includes the n-color simple exclusion process, which is a standard simple exclusion process 
but each particle is colored with one of n possible colors. Exchanges only occur between 
occupied and empty sites, so that two sites with different colors cannot exchange. Our 
arguments do apply to this case too. This can be viewed as a multi-species system and the 
effect of long range interactions is of interest in its own right and we discuss it in Section 
5 of this paper, limiting ourselves to the n = 2 case, corresponding to a binary alloy in 
which exchanges take place only through vacancies. 

The paper is organized as follows: in section 2 we give the precise definition of our 
model and we state the main results. In section 3 we recall some important results of 
[VY] that we need for our proof. Section 4 contains the proofs which are based on the 
control of the Radon-Nykodim derivative of the full process with respect to the reference 
process and on the [GPV]-[V] method. Finally in section 5 we consider the example, that 
we just mentioned, of a dynamics with two conservation laws, for which we prove a result 
analogous to the one in section 2. Moreover, in this case we get a stronger convergence 
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result, since uniqueness of the weak solution for the limit equation does hold, under suitable 
hypotheses. 

2. The model and the main result. 

We work in the discrete torus of dimension d and diameter 7 _1 , 7 > 0, that we denote 
by A T . Associated to A 7 there is a natural notion of A* , the set of (non directed) bonds, i.e. 

the couples of nearest neighbor points of A 7 . The configuration space on A 7 is {0, 1} A ^ = 
0a 7 - All these spaces are endowed with the discrete topology: when later we will deal 
with {0, 1} A , A countable, we will use the product topology. 

2.1 The Gibbsian reference measure. Let F be a local isotropic function from Z d 
to IR. By isotropic we mean that F(9rj) = F{rj) for every reflection 9 with respect of an 
axis as well any lattice rotation 9. For any x G 7L d set F x (rj) = F{r x rj), where t x is the 
translation of x G ^ d in Z d : 

TxViv) = V(x + y), for every y G % d . 

And r will also be the translation operator acting on functions F of the configuration: 
t x F( V ) = F(t xV ). 

Let us consider the formal Hamiltonian 

H{r } ) = Y J Fx{rf). (2.1) 

X 

It is well know that if the inverse temperature k is sufficiently small, given any chemical 
potential A, there exists a unique extremal Gibbs measure p p , p G [0, 1], which satisfies p = 

E /iP [rj(x)] for every x G 7L d . We will always work in this uniqueness regime: this guaranties 
also that p p is translation invariant and that it has some mixing properties. However in 
some proofs in [VY] the authors require a stronger mixing condition (exponential mixing 
condition, [Assumption A, VY]), given in terms of finite volume grand-canonical measures. 
Since we are using this assumption only indirectly we do not give it explicitly: we just 
stress that our results are proven only if k is smaller then a certain kq > 0, which depends 
on the dimension d and on the interaction F. In some cases it can be shown that kq = k-c, 
the inverse of the critical temperature of the reference system. 

We will deal with the F-interaction also in the case of A 7 : for 7 sufficiently small we 

can keep the very same infinite volume definitions by lifting A 7 to the whole of 7L d in the 
natural way. This way the notion of translation t x is unchanged. It is however often more 
natural to look at the periodic finite volume case, and that's what we will do. 

The same applies to every local function on Z d, which will be viewed as a function on 
0a 7 without notational changes. 

2.2 The equilibrium measure for the full system. We consider a the probability 
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measure on Oa denned by 



= z ( J iKjX) ex p { -« E F * - § E J 7 (*« + A E rt*) } ' (2-2) 

13 > 0, 7 > 0, A G R, 

where J 1 {x,y) = ^ d J(^(x - y)), J G C 2 (T d ) IR (T d is the d-dimensional torus of 
diameter 1) such that 

J(r) = J(-r), [ J(r)dr=l, (2.3) 

and Z^((3, k, A) is the normalization factor. 

2.3 The dynamics. This will be introduced, for the reference system, both in the case 
A = A 7 and in the case A = Z d . Let $ G C 2 (IR; IR+) be of the form 1 

*(E) = e X p(-E/2)<l>(E), 

(2.4) 

< />:IR-^IR + , <j>{E) = <j>{-E) for every E G H and 0(0) = 1. 
For b = (x, ?/) G A* we define 

c°(r ? ) = $( K A b if(r 7 )), (2.5) 
where A b H(r]) = H(r] b ) - H{rj) and 

rj(x), iiz = y, 
ri b (z) = { V (y), ifz = x, (2.6) 
rj{z), otherwise. 

In the infinite volume case H is not well defined, but A b H is taken, by definition, to be 
equal to lim^^oo A b Y1\ x \<r Fx(v)- To make the notation a bit lighter, if b = (x, y) appears 

as a subscript, we will often drop the brackets. Given / : Oa — > IR for b = (x, y) we set 

C° b = c° b (v)[f(v b )-m}- (2-7) 
A Markov pregenerator is then defined by setting 

C°m = J2 £V(V), (2-8) 
beA* 



1 This is equivalent to the more customary [p. 163, Sp] detailed balance condition <&(E) 
$(-E)exp(-E), E G IR. 



where / is assumed to be local in the case A = If A = A 7 , C° is actually a generator 
and it is easy to construct a unique process in Skorohod space {ilt} t en+ G D(IR + ;0a 7 ) 
associated to it, once an initial condition is given. The law of this process will be denoted 
by P 7 ' or P 7 ^ if there is the need to stress the initial condition p 1 . It is immediate to 

verify that, by (2.4), C° viewed as an operator in L 2 (p^ K,x ) is self-adjoint for every A G H. 
In the infinite volume setting the construction of the process is more delicate. We refer to 
[Li] for this construction: there the process {f]t\teR + G -D(M + ;0 Z d) associated to C° is 
constructed. We remark also that, for every p G [0, 1], C° can be extended to a self-adjoint 
operator on L 2 (p p ). The law of the process {?7t}t 6 R+ will be denoted by P°. We will 
sometimes stress the chosen initial condition, say p G Pi (Ha), by writing P°. Here we 
used Vi(-) to denote the probability measures on •. 

The full dynamics is considered only in the case A = A 7 . We define 



K 



x,y 



(2.9) 



and 



£y/fa) = £ <£fa) [f(v b )-f(v)} 

beA* 



(2.10) 



As before, associated to the finite dimensional operator C 1 there is a process with trajec- 
tories in D(IR + ;0a 7 )- £ 7 is self-adjoint in L 2 (p^' K ' x ) for every A G IR+. The law of this 
process will be denoted by P 7 = P^ 7 , where p 1 G V(£l\^) is the initial condition. 

2.4 The main result. The compressibility x f° r the system is defined as 



X (p)= £ cov^^O),^)), 
xez d 



(2.11) 



in terms of the local interaction alone. The corresponding diffusion matrix D can be 
expressed via a variational formula. It is the p dependent symmetric matrix defined by 



(v,Dv) M d 



inflE^ 

2X(P) a 



£4(7?) j vMa) -77(0)) - £ A (0 , ei) g(T xV ) 



i=l 



(2.12) 

for every v G R d . In (2.12) denotes the unit vector in the % direction, (•, -) U d is the 
scalar product in IR d and the infimum is taken over all local functions g. We will take the 
freedom of using both the notation Vi and the notation v ei . Moreover for e a unit vector, 
v e = ±v ei if e = ±e;. 

A number of facts are known about D: first of all it is a continuous function of p, cf. 
[VY] , and there exists a constant c such that in the sense of matrices 



< D(p) < cl, 



(2.13) 
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where I is the dx d identity matrix. While the upper bound is an immediate consequence 
of (2.12), the lower bound is much more subtle and it has been established in [SY]. In [VY, 
Lemma 8.3] it is shown moreover that in our isotropic set up, D(p) is a multiple of /. We 
will keep the notation D(p) also to denote the scalar proportionality factor between / and 
the matrix D(p). 

We are going to consider initial particle configurations associated to a density profile 
Po : T — >• [0, 1], in the following sense: if we define the the empirical density field 

^(t, x) = 1 d Y. 6 ( x - w)^-n(y), (2-14) 

yeA 

we require that for any smooth test function G from T d to IR and 5 > 



lim P 7 

7^0 ^ 



z/ 7 (x, 0)G(x)dx — / po(x)G(x)dx 

'T d JT d 



>5) =0. (2.15) 



Denote by Q 7 the law of the process {^ 7 (t)}te[o,T] on the space D([0,T\,M), induced 

by P^ 7 , where M. = Ai(T d ) is the space of nonnegative measures on the torus with 
total mass bounded by 1: this space is endowed with the weak topology weakened by the 
continuous functions. We consider also the subspace M\ = M\(T d ) of M. consisting of 
absolutely continuous measures with density bounded above by one. Our notation for the 
gradient in <i-dimensional Euclidean space is d. 

Theorem 2.1 Consider an initial datum satisfying (2.15). Then, the sequence of probabil- 
ity measures Q 7 is tight and all its limit points Q are concentrated on absolutely continuous 

paths whose densities p e C°([0, T]; Mi(T d )) n L 2 ([0, T], #i(T d )) are weak solutions of 

the equation 

d t p = d{D(p) [dp-0 X (p){dJ*p)]} , 

(2.16) 

P(0,-) = Po(-). 

Moreover if the diffusion coefficient D is Lipschitz continuous, then Q 7 converges and the 
limit point Q is the unique weak solution of (2.16). 

Throughout the text <E>, J, k and (3 will be considered fixed. We will be often interested 
in getting estimates which are uniform on the configuration n or on the history of the 
process {rjt}t>o- We therefore introduce the notation o u (-) and O u (-) in the standard 
sense of o(-) and O(-) but uniformly with respect to the configuration or to the history of 
the process. 

3. The fluctuation— dissipation equation. 

In this section we recall the fundamental result proven in [VY]: the approximate de- 
composition of the current in a gradient term and a fluctuation term. 
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Let us recall the definition of current in the general J-dependent case. The current is 
defined for every x, every unit vector e and every rj as 

:+ e (v) = cl jX+e (v)[v(x) ~ V(x + e)]. (3.1) 



w 
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It has the property that 

d 

C^ix) = -J2 [ w 2,x+ ej (v) - w2- ej , x (v) • (3-2) 
i=i 

The analogous definitions in the case of the unperturbed dynamics generated by C° are 
just obtained by setting 7 = 0. 

We follow very closely [VY]. For ( e n zd or ( : Z d — > R and £ G IR+ we set 

V ; y-\y-x\<l 

We also set A\ = {x G Z d : |x| < £}. We define the a-algebra J- X ,a generated by 
{A\ s r](x)} U {rj(y) : \y — x\ > s} and the space Q of local functions h : Z d — > 1R with 
the property lE^f/il^] = for some s. To any h G Q and any 77 G Z d we associate an 
element (still denoted by h) of Q Z d by setting h(x) = T x h(rj). We call Vtg the subset of 
Z d obtained from with this procedure. We then introduce the finite volume variance 

V t (h,p,t) =£ d (Av ei h,(-£° A ,)- l Av ei h) flA/iP ^ (3.4) 

where £\ = £ — VI, £° = Xl&eA* ^6 f° r an y finite set A C Z d and PA' e ,p,£ is the canonical 

Gibbs measure with interaction F on A^ with boundary condition £ G Z d and density 
p G [0,1]. Two observations are in order for canonical measures. The first is that we 
view PA' e ,p,£ as an element of the probability measures over Z d: the extension is made in 

the natural way obtaining a measure concentrated on {rj : rj(y) = £(y) for every y G A' £ c }. 
Second, if p is not an integer multiple of 1/|A^|, the density of the canonical Gibbs measure 
is taken to be max{/c/|A^| : k G ^, k < p\A' £ \}. We observe then that if x G , h(x) 
depends only on rj(y), y G A' e , for sufficiently large £. Moreover Vg(h, p, £) is an Tt- 
measurable function of £ G O. Given pP ', the infinite volume Gibbs measure with density 
p, we define V : x [0, 1] -> IR + by 

y(/i,p) = limsup / (^,P, CM/AO)- (3-5) 

Finally we extend this definition to every local function h by setting h}-(x) = h(x) — 
E MA r"' f [h{x)\F x , k \ and taking the limit 

limsup V(hk, p)- (3-6) 

fc^oo 

Below we use the notation 

V ej r)(x) = rj(x + ej) - rj(x). (3.7) 
We are now ready to state the Fluctuation-Dissipation Theorem [VY, Theorem 3.4]. 
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Theorem 3.2 The symmetric, in fact diagonal, density dependent matrix D defined in 
(2.12) coincides with the following Green-Kubo formula: 



X (p)D iJ (p) = 8 id hE^ 



'(0,e,-) 



/•OO 



w o, ei ex P t x w% dt, (3. 



where x is the compressibility defined in (2.11) and Sij is the Kronecker delta. Then for 
any a G JR d 



/ a 



P =0. 



(3.9) 



Moreover, for any 6 > t/iere exists o 5 : [0, 1] x Z d — > IR d , Q 5 (p, •) G £ d /or every p and 
g s (-,i]) is smooth for every n, such that 



sup V f J2<*j \wleM + ( D b)Vv(P))i + £°9 
pe[o,i] T L 



p <5. 



(3.10) 



We observe that, by polarization, from the bilinear functional V(-,p) we can define a 
scalar product (covariance) that will be denoted by ((f,g))(p), defined for / and g local 

function on {0, 1} Z . This covariance is carefully analyzed in Section 8 of [VY]. We collect 
here two properties that will be crucial for us. First, formula (8.7) in [VY] tells us that 

(Ke,^o°e,»(p) = \&ij& P [co, e ,(r/)(V e ,r7(0)) 2 ] . (3.11) 

Moreover formula (8.13) in [VY]) tells us that 

(K°e,-V,#)))W = ^x(p). (3.12) 



4. Proof of Theorem 2.1 

4-1. Preliminary lemmas. We will repeatedly need the expansion (in powers of 7) of the 
jump rates: we take advantage of the fact that an exchange of two particles changes the 
long-range energy of O u (l)- We will use the following notation for discrete convolution 

(f*v)(x)=Y E /(7(*-*))»7(*), (4-1) 
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where / is a function from T to IR and x G A 7 . 

Lemma 4.3 For every $ and J, there exists C such that for every b G A 7 *, every r\ G fl\ 
and every 7 G (0, 1) 

A 6 ( J 7 (x,y)rj(x)ri(y) 

and 

1^(77) - c° b (n) - 7/3$' («A 6 tf (77)) [ V (x + e) - t/(x)] e • (9 J * 77) (x) \ < C7 2 , (4.3) 
where b = (x, x + e) . 

Proof. The proof follows immediately from the expansion in Taylor series of $ (2.4). □ 

A first application of Lemma 4.3 is in proving that the concept of superexponentially 
small event coincides for perturbed and unperturbed processes. As usual, (3 and k are 
fixed, but recall that we denote by P 7 ' the law of the process in the box A 7 with (3 = 0, 
that is the reference process in the periodic box. 

Lemma 4.4 Let T be a fixed positive number, A = A 7 and let £> 7jq be a bounded func- 
tional of the process {r]s}se[o,T'y- 2 ] which depends on 7 > and on a vector parameter q. 

Furthermore let fi = fj,® ,K ' (recall (2.2)) and let fi' be any probability measure on 0a 7 - If 
for every p G IR 

lim lim sup 7 d log E 7 ' (exp {p^~ d B 1 ^ ) = 0, (4.4) 

then the same is true with E 7 ' replaced by E^, . Analogously if {i? 7i q} 7i q is a family of 
J-'r-y- 2 measurable events which are superexponentially small under P^'°, that is 

lim lim sup Y log P^'° (£ 7 , q ) = -00, (4.5) 
then they are superexponentially small also under P^, . 

Proof. First we observe that, since for some c = c(F, k, A) G R we have that sup M , n (d/z'/d/i) (77) < 

e c|A 7 |^ ^ - g su ffi c i en t to prove the statements in the case fs! = fi. We therefore omit the 
subscript fi in this proof. 

For b = (x,y) and t > 0, let iV&(t) denote the number of jumps between a; and j/ in 
the time span [0,t]. The Radon-Nikodym derivative of P 7 with respect to P 7 ' , both 
restricted to J-'t-y- 2 wm be denoted by M t and it is given by 

M t = expf-jf* 7 ^KM-cgfa^da + jT £ log di\fc( a )) , (4.6) 



<C 7 , 



(4.2) 
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c.f. [Prop. 2.6, App. 1, KL]. With respect to the measure P 7 ' and the filtration {jF t } t > , 
the process {M t } t >o is a martingale. Therefore, for p > 1, {Mf }t> is a submartingale. If 
we define 

A t = (4.7) 



f 1 pM£j2l c UVs)-c° b (Vs)}ds+ F M^^n 
Jo Jo 



plog' - i 



ds, (4.8) 



we have that Mt = Mf — At is a martingale with Mo = 1. We now expand the expression 
in the right-hand side of (4.8) by taking advantage of the fact that (cZ(r} s )/c®(r} s )) — 1 = 
0^(7). Precisely, by Lemma 4.3, 



exp ( p log 



c° b (Vs) 



1 = p 



C° b (Vs) 



1 + 



C° b (Vs) 



1 +O u ( 7 J ), (4.9) 



where we have assumed p bounded, say p < 2, for the last term. Therefore we obtain that 
there exists C such that 



< C{p{p - l)7" d+2 + 7" d+3 ) / Mf ds, 

Jo 



(7 3 ) 



ds 



(4.10) 



where in the last step we have used the fact that c^rj) is uniformly bounded as well as 
the positivity of M t . By taking the expectation of this last expression, recalling that 
E 7 '°[Mt] = 1, by Gronwall's Lemma we obtain 



E 7 ' [Mf] < exp (Ct (p(p - l) 7 " d + 7" d+1 )) . 



(4.11) 



This suffices for our purposes since, by applying Holder inequality, we obtain that for every 
p > 1 and q = p/(p — 1) 

lim sup lim sup 7^ log E 7 (exp(p7 _d £> 7)q ) 

< lim sup — logE 7 ' (M£) + lim lim sup — logE 7 ' (ex.p(q~ f - d B 1 q ) 

7^0 v 1 7^0 q 

< Ct(p-l). 

(4.12) 

By letting p \ 1 the first statement is proven. The proof of the second statement runs in 
the same way. □ 

If h is a local function, we define h{p) = lE /iP [h] for p E [0, 1]. Here is the first application 
of Lemma 4.4. 
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Lemma 4.5 {Replacement Lemma). Let h be a local function and 



B bl -i{rj) 



1 



(2&7- 1 + l) d 



[h(T v v) ~ h (Av [fry -i]7/(0)) 



y\<b~i~ 



(4.13) 



for b > 0. Then, for any 5 > 



lim sup lim P 7 



Jo 



= 0. 



(4.14) 



Proof. This goes through the by now classical one block and two blocks estimates. These 
can be found in [VY] (Lemma 5.2 and Theorem 6.2) for the unperturbed process and these 
estimates are superexponential. The extension is therefore just Lemma 4.4. □ 

We conclude this subsection with a computation that is very relevant for us to identify 
the limit equation. 

Lemma 4.6 . For any bounded local function f : Z d — > 1R we have that 



[^( K A , e H)V eV (0)A ,ef(v)] = ~\^ P m^A 0je H)V eV (0)A , e f(v)] • (4-15) 



Proof. By differentiating both sides of (2.4), we reduce (4.15) to proving that 
E » P [ eX p (-«A 0l e#/2) 0'(KAo, e #) (77(e) - 7/(0)) A , e f(r,)] = 0. 



(4.16) 



Let us now approximate \x p with a sequence of finite volume grand-canonical Gibbs mea- 
sures on 0a 7 . The result follows because 

= £ e-* A °.^)/y(«Ao,e#(r7)) (77(e) -^(0)) Ao.e/^e-^e^^W 



£ (77(e) -7/(0)) Ao, e /(7 ? )e- K ^W + ^ ' e )]/V / («A , e if(7 ? ))e A S^(-) 



(4.17) 



£ (77(e)- 77(0)) Ao, e /(7 ? )e-^ ' e )+^/ 2 ^(-/,Ao, e ^(7 ? ))e A E,^)^ 



where in the last step we used the fact that (77(e) — 77(0)) Aq jB f(r]) is invariant under the 
transformation 77 — > n 0,e . Recall now that 4>' is odd and the proof of (4.15) is complete. 
□ 
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4.2. Tightness and energy estimate. 

We go now to the set up of Theorem 2.1. Recall that Q 7 is the law of the empirical 
process {u 1 (t)}t£[o,T], cf. (2.14), on the Skorohod space D([0,T], M). 

Lemma 4.7 The sequence {Q 7 } 7 >o is tight and every limit point Q is concentrated on 
absolutely continuous paths whose densities belong to 
C° ( [0, T] ; M i (T d ) ) n L 2 ( [0, T] , H x (T d ) ) . 

Proof. This is standard. One way to prove tightness in D([Q,T],JA) is to repeat the 
proof in [§4, VY]: Lemma 4.1 in [VY] depends only on the uniform boundedness of the 
jump rates and Lemma 4.2, still in [VY], is easily upgraded to our situation via Lemma 
4.4. The limit points actually lie in a smaller space: by the exclusion rule it is immediate 
to verify that the limit is in Ai± for every t and, since every jump produces a discontinuity 
O u (l d ), we can substitute D with C°. The energy estimate, that is the existence of a 
constant C such that any limit point Q satisfies 



(Vp(s, r))drds 



< C, (4.18) 



requires a more sophisticated argument. Once again most of the work has been done in 
[VY], Section 5: it is sufficient to upgrade (5.25) and (5.28)-(5.30) in [VY] to our situation. 
One way would be to get a volume order estimate on the relative entropy of P 7 with respect 
to P°, both measures restricted to [0, T7 2 ], but we choose to stick to LP estimates on 
the Radon-Nicodym derivative Mt of the process. We will be as close as possible to the 
notations of Lemma 4.4 and, like in its proof, we will omit the dependence on the initial 
condition: the only difference is that we will not need any extra parameter q and we will 
therefore omit the subscript q from the notation. We choose, with G G C 1 (R + ,T d ;R) 
and C\ > 0, 

£ 7 , q = B 1 = Y' 1 [ ^G(t,x 7 )V e ^ 7 - 2 (x)dt-^ 7 d I ^T|G(WiV)| 2 dt (4-19) 
J o x 1 J o x 

In the proof of Lemma 3.2 in [VY] it is shown that for some C\ (4.4) holds for p = 1, that 



is 



lim 7 d logE 7 ' [exp (7"%)] = 0. (4.20) 

Using (4.20) we will show that there exists C2 such that 

limsupE 7 [B 7 ] < C 2 . (4.21) 

We refer the reader to [Section 5,VY] for a proof of the fact that (4.21), with the choice 
(4.19), implies (4.18) with C = max(Ci, C2). Here we prove (4.21); it is just an application 
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of Jensen's inequality and Cauchy-Schwarz: 



E 7 [B 7 ] < 2 7 d logE 7 ' 



M^exp 



7" 



(4.22) 



< 7 d logE 7 ' [M|] + 7 d logE 7 ' [exp (7"%)] , 

and, recalling (4.20), (4.21) follows from (4.12) with p = 2. □ 

4.3. The (partial) identification of the limit. By the definition of the process, for 
every G E (^(T^H) we can write 



Jo 



(4.23) 



in which = {M^(t)} reU + is a P 7 -martingale with respect to the filtration associated 

to {r]t}teu+- The quadratic variation of is easily computed and estimated: if we set 

Xj(v) = l d J2 X G{ix)r](x) it is immediate to see that there exists a constant c = c{d,G) 
such that for every 7 > 

ft-/' 2 

(M°M?)(t) = / (C 1 X 2 -2X 1 C 1 X 1 )(i ls )ds<c 1 d , (4.24) 
Jo 

and thus, by Doob's inequality, for every T > and every 5 > 

limP 7 ( sup \M^(t)\ > 5 ) = 0. 



7^0 



t€[0,T] 



(4.25) 



We now split the current: 

w t = [ w °] + [ w -r - ( 4 _ 26 ) 

The non gradient difficulties come from the first term. Below we use the convention that 
J2 e is the sum over the unit vectors {ej}j = \^,...,d- 

We first have a fast look at the easy second term. By Lemma 4.3, for every x and every 
unit vector e 

[w^ - w°] x , x+e (rj) = -/? 7 $' (A x , x+e H( V )) [V eV ] 2 (x)e ■ (dJ* V )(x) + O m ( 7 2 ). (4.27) 
If we define 

Kb,-, (R} s > ) = j 7 d E [G(7(« + e)) " G(7*)] [< I+ e(lt) - <** 
+/?f V +2 E (g e G) (7*) (d e J * Vt) (x) . I )d E Mwfc)d* 



(4.28) 
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where h e (rj) = 3>'(Ao je H(r}))(r}(e) — ??(0)) 2 , and therefore 



h e (p) = [$' (KA , e H( V )) (77(e) - V (0)) 2 } . 
Note that, by applying Lemma 4.6, with / = 1, we can immediately rewrite 



(4.29) 



h e (p) = --W P [c° 0t M(v(e)-v(0)) 2 } 
By the smoothness of G and J it is immediate to obtain that 

limlimsup sup TZ bn [{r] s } s>0 ) = 0. 
b-o 7 ^ { Vs } 3>0 V - / 



(4.30) 



(4.31) 



This, together with Lemma 4.5 applied to the spatial average of h, immediately implies 
that 



/ 7 d ^[G(7(£ + e)) -G(7^)] H, x+e {m) -<x+eM dt- 



-/T d 



dG(r)(8J ★ z/ 7 )(t, r)/i (l ( _ 6j6)d * z/ 7 )(t, r)) drdt 



(4.32) 



tends to zero in P 7 -probability as 7 — > and then b ^ 0. 

Now we turn to the non gradient term. Given a, 6, 5 and <5i > 0, we define 



Ba,b,5 = <i SUp 

te[o,T] 



7 



d+l 



7" 2 t 



^^d e G( 7 *)< x+e (^)dt+ 



1 / ( 26 ) _1 l-Vsix + 6 7 _1 e) - Vs(x - &7 _1 e)] £>e,e' (Av [a7 -i]7/ fl (x)) a e G( 7 x) + 

7 <m r T ^a e G( 7 x)£V 5 (^^) 



7 



><*i 



(4.33) 

The following result follows from Theorem 3.3 in [VY] and Lemma 4.4, with the only 
observation that in the statement in [VY] the term C°g s does not appear: in their case it 
is irrelevant (see Section 3 of [VY]). 



Lemma 4.8 For every 8\ > 



lim lim sup lim sup lim sup P 7 (B a ,b,s) = 0. 

6—0 <5— 7—0 



(4.34) 
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We now observe that in (4.33) we can replace C°g s with (£° — C 7 )g s and Lemma 4.8 
still holds since 



7 



d+l 



x,e 



1 



7 d d * G M (9 6 e (^-zt) - 9 S e (r xVo )) 

x,e 



^d e G(<yx)C 1 g s e (r x r]t) 

+ 7 M 7 (T) + £ 7 (T). 



(4.35) 



where {M J (t)}t>o with respect to the natural filtration is a P 7 -martingale and E^(T) is 
the error we made by considering g s independent of {Avg s rj(x)} x , is the range of g 5 . The 
first term in the second line of (4.35) is u (i), the martingale M 7 has a quadratic variation 
0(7 d ) and sup t€ [ T ] |.E 7 (£)| tends to zero in probability as 7 — > and 5^0. The last 
claim follows from the argument in [VY], formulas (3.26)-(3.27), and Lemma 4.4. 
We are therefore left with evaluating 



7 / i d d eGM [C° - £ 7 ] g s e (r xVt ) dt. 



(4.36) 



By using Lemma 4.3 to express C° — £ 7 , recalling that g 5 e is a local function and repeating 
exactly the same steps, i.e. smearing and using the replacement lemma, as in (4.28)-(4.32), 
we obtain that, up to a negligible error term, this term can be replaced by 

pY" I I d e G(r)6l e ,( X[ -6,6]**^(t,r))d el J*v'i(t,r)drdt, (4.37) 



where 



E + e') - rj(z)) $' (KA z , z+e ,H ( v )) \g s e (v z ' z+e ') - g S e (v) 



. (4.38) 



Let us define 9{p) = limsup^g 9 s (p) for every p e [0, 1]. By Lemma 4.7 and collecting 
the results (4.23), (4.25), (4.31) and Lemma 4.8, under the hypothesis that 9 s converges 
uniformly as 5 — > 0, we obtain that every limit of the sequence {^ 7 } 7 >o concentrates on 
trajectories p which solve weakly the PDE 



dtp = J2 di { D id(p) d iP ~ P (kjhjip) + 9i,j(p)^ (djJ * p)} 



■1 -j 



(4.39). 



P(0,-) = Po(-). 
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We are therefore left with showing that 

Lemma 4.9 The sequence 9 s converges uniformly and for every p E [0, 1], every i and 
every j = 1, 2, . . . , d we have that 



$i,jhj(p) + d i,j(p) = x(p)Di,j(p). 



(4.40) 



Proof. We divide this algebraic computations in steps 

Step 1. We start by rewriting 9 in a more convenient form, using Lemma 4.6 with / = gf . 
We obtain that 



J2 c z,z+ ej {v) (v(z + ej) - rj(z)) A Z!Z+ej g s e .(ri) 



(4.41) 



Step 2. Reversibility and summation by parts on (4.41) imply 



0tj(p)=^ 



= -Ke 3 .>Se t >0(p). (4-42) 



where in the last step we used the notation 



(/,<7>o(p) = J>^ \f(v)g(r x v)] , 



(4.43) 



for / and g local functions. It is not too difficult to show [Section 8, VY] that 

{f,g)o = -((f,£og)h (4-44) 
and we recall that the covariance ((■, •)) is defined right after Theorem 3.2. Therefore 

0l j = («e j ,£o9 S e i )). (4.45) 



Step 3. We now use the decomposition induced by (3.9) to express Cog 5 e . obtaining that 
6f j is equal to 

= (K^OIW - (K^WJOj^J + ^i), (4.46) 

where lim^o sup p \R(p, S) \ = 0. 

Completing the proof of (4.40) is now just a matter of applying (3.11), note the cancel- 
lation with the term h, c.f. (4.30), and (3.12) and the proof is complete. □ 
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Lemma 4.9 leads us to the weak formulation of the PDE (2.16), that is that every limit 
measure Q is concentrated on trajectories p G C°([0,T];Mi(T d )) n L 2 ([0, T], Ui(T d )) 
which solve 

/ G(r)p(T,r)dr- I G(r)p(0,r)dr = 

JT d JT d 

[ [ dG(r) (D(p(t, r)) X (p(t, r))dJ * p - D(p)dp(t, r)) drdt, 

JO JT d 

(4.47) 

for every G G C^T ') and every T > 0. Therefore the proof of Theorem 2.1 is complete 
modulo discussing the uniqueness issue. This is considered in the next subsection. 

4.4. Uniqueness. A proof of uniqueness is available if D is a uniformly Lipschitz 
continuous function. This follows from the control of H~ x norm. In [GL2] this result is 
proven under an additional condition on the time derivative, in the sense of distributions, 
of p. This condition can be removed if one replaces the kernel of the H~ x norm with a 
smoothened kernel, as it is done in [Appendix A, KL]. Note that the weak formulation 
(4.47) is in this set up equivalent to the formulation that we get by multiplying both sides 
of (2.16) by G G C 1 (IR + x T d ; IR) and formally integrating by parts. 

5. Multispecies systems 

The scheme of proof of Section 3 and Section 4 can be applied to several other systems. 
Of particular interest for applications is the case of systems with several types of particles. 
One (apparently) simple system in this class is the n-color exclusion process: take a simple 
exclusion process and distinguish the particles by painting them with n colors (n is kept 
fix). The hydrodynamics of this system has been done in [Q]: here we would like to consider 
the case in which the particles interact via long range potentials that distinguish between 
colors. 

From the technical viewpoint this case has essentially been already considered in [QRV] , 
where a Large Deviation principle for a n-color exclusion system is proven as a step to 
prove a process Large Deviation for the simple exclusion. The lower bound of the Large 
Deviations depends on the standard change of measure argument which entails studying 
the n-color system driven by a weak external force. This is also our case: the weak driving 
force is however configuration dependent, but since it depends on the configuration only 
via an empirical average, the changes with respect to [QRV] are minimal. We will therefore 
simply state the result and make some observations on the proof. 

To simplify the notation and the statement of the result we restrict our attention to the 
case of two species (A and B) and to the case in which A and B interact with each other, 
but A does not interact with any A and the same for B particles. 

5.1 The A— B model. As in the previous sections A 7 will denote the lattice torus with 
diameter [7 -1 ]. We are looking at a <i-dimensional system of A and B particles evolving 
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via a Kawasaki dynamics with Kac Hamiltonian 

H i(v) = l E Mx,y)v A (x)v B (y) (5.i) 

where as usual 7 > and ry" 4 and ?y B are elements of {0, 1} A ^ with the hard core restriction 
that there can be at most one particle per site 

V A + V B <1. (5.2) 

J 7 (x, y) = "f d J{l{x — y)), where J is a smooth function from the <i-dimensional unit torus 
to IR. The particle configuration can be alternatively described by 

77 G {0,^4,5}^, (5-3) 

and we will identify 77 with (r? A , 1] B ). The dynamics which conserves both A and B particles 
(or interchanging with) are specified by A and B particles hopping to nearest neighbor 

empty sites at a rate Dc^ ,AB (rj) but no direct exchanges between A and B particles are 
permitted. This models the physical situation of polymers in a fluid considered in [WD]. 
The generator of the dynamics is 

L 7 /fa) = (5.4) 

d £ c 2 ,AB (v) [(v A (x) - v B (y)) 2 + (v B (x) - v A (y)) 2 } if(v x ' y ) - f(v)} , 

b=(x,y)EA* 

where D > 0, / is a real valued function defined in {0, A, B} A f and i] x ' y is 1] with the 

occupation of the sites x and y exchanged. As before, the rates cl' AB are such that the 
dynamics satisfies the detailed balance condition with respect to the Gibbs measure with 
Hamiltonian (5.1) at inverse temperature (3 > 0: 

c (^f fa) = $ ^ l H ^ y ) - H y ( V )]) , (5.5) 

$ as in (2.4). To denote the process we keep the same notation of the previous sections. 

5.2 The hydrodynamics of the A— B model. For the hydrodynamic limit we look at 
the empirical measure 

"2(*> x) = 7 d E 6 ( x ~ TVH-'M ( 5 - 6 ) 
yeA-y 

for a G {A, B}, x G T d and r] a (x, s) specifies the presence or absence of an a-particle on 
site x at time s. 

For the initial datum we assume ^ (0, •) to be close to p Q (-) in the sense of (2.15). In 
this case we have to make further assumptions on p Q , precisely that there exists 5 > such 
that for every x and every a 

5<p a (x)<l-5, (5.7) 
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and that p a is differentiable with bounded derivatives 

sup |<9p a (;r)| < oo. 



(5.8) 



x,a 



For the moment let us set J = 0. In [Q] it is shown that converges weakly in 
probability, that is in the sense of (2.15) for every t > 0, to p a : IR + x T d — > [0, 1] and the 
couple (p A 7 p B ) solves 



dt 



= d- 



B A 

p s p 

„S . A 



^{D-D s ) 



d 



-A 



(5.9) 



£-(D-D a ) ^D S + ^D 

p v */ p p 

with initial condition (p A (0, -),p B (0, •)) = (p A (-), p B (-)). In (5.9) p = p A + p B and D s = 
D s (p) is the self diffusion coefficient. The expression for the diffusion matrix D in (5.9) 
can be derived from elementary considerations on the microscopic system from which it is 
derived [LS]. We observe that, as expected, the evolution equation for p is simply 

dtp = DAp. (5.10) 

This follows from the observation that if we ignore the distinction between A and B 
particles then, in the absence of interactions, we are just dealing with the one component 
simple symmetric exclusion process [Sp]. 
We can rewrite the system (5.9) as 



d t p = V- [VVp] = V 



MV 



in which p = (p A , p )*, and 



T I A B\ 1 f / A B 

Fo{P ,P ) = ~p J s{p ,p 



Sp 



)dx, 



with 

and 
Explicitely 

and 



s(p A , p B ) = - [p A log / + p B \ogp B + (l-p A - p B ) log(l - p A - p B )] 

M = -/3P(Hess(s))- 1 . 



Hess(s)- 1 = 



i_(-P A (1-P A ) 



P A P B 



P A P B 
-p B (l-p B ) 



M = /3 



Dq 2^1 + D (P A ) 2 (^) 
p p 

^B„A n „-> A n B 



s p 

B\2, 



D 



p a p b (i-p) 



(5.11) 

(5.12) 

(5.13) 
(5.14) 

(5.15) 
(5.16) 



We now claim that, in the case in which J is not zero, the limit equation (5.11) has to 
be changed by replacing JF with 



F=F +^J J J(x-x')p A (x)p B (x')dxdx'. 



(5.17) 
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Proposition 5.10 Set d > 3. For every t > the empirical field (z/J(t, •), ^(^ ')) 
converges weakly in probability, i.e. in the sense of (2.15), to the unique weak solution of 

d t p = V • [Vdp + Md J * p] . (5.18) 



The restriction to d > 3 is due to the fact that only under this restriction D s is known 
to be a Lipschitz function and uniqueness of the weak solution holds. 
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